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ABSTRACT 

This thesis considers the problem of providing attitude 
control for a spacecraft engaged in an extended mission. As a 
basis for the choice of a suitable attitude control system the 
following requirements are applied. 

Maximum reliability 
Minimum ejection of mass 
Minimum average power 
Minimum system weight 
Minimum peak power 

An interplanetary mission of 400 days duration is adopt- 
ed as a general guide for the problem, but most of the equations 
and comparisons are presented in parametric form. Extended 
missions imply that a momentum exchange type attitude control 
system be used to minimize ejection of fuel mass, and the thesis 
primarily considers only systems of this type. The thesis derives 
the equations of motion for a spacecraft equipped with eighteen 
different control systems. The control system chosen to best 
satisfy the five design requirements is a system consisting of 
four gyro-type controllers arranged in two pairs with each pair 
operating back-to-back to minimize control cross coupling 
torques. One pair of controllers provide roll torques, the other 
pair provides pitch torques, and all four controllers contribute 
yaw torques. 

The stability and control analysis considers operation 
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contribute yaw torques. 

The stability and control analysis considers operation 
of the spacecraft in three modes. 

Zero Input Mode 
Rate Control Mode 
Position Control Mode 

Each of the modes are evaluated for roll motion by assuming 
negligible interaxial coupling, and the analysis includes opera- 
tion of the controller gimbal angles to large angles. 

Thesis Supervisor: Wallace E. Vander Velde 

Title: Associate Professor of 

Aeronautics and Astronautics 
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CHAPTER 6 



CLOSING THE ATTITUDE CONTROL LOOP 



6. 1 Introduction 

Previous chapters have defined the physical character- 
istics of the spacecraft, have chosen the four gyro (12-34-12 34) 
system as being the most attractive control system, and have 
defined probable disturbances on the spacecraft. Before the 
mathematical models representing the spacecraft and control 
system can be combined into a single loop it is necessary to 
specify an attitude sensor and some type of torque motor which 
will receive signal information from the sensor and provide the 
proper angular rates and positions to the spin axes of the four 
controllers. The sequence of this chapter follows an order that 
first determines the dynamics of the individual controllers by 
methods similar to those of Chapter 2 which considered the com- 
plete spacecraft. Next, the motors and drive systems which 
drive the controllers are defined. Then, the sensor and related 
amplifiers are assumed together with the provision for tandem 
compensation, if required. Finally, all of the above components 
are combined into a single loop. 

6.2 Response of the Controllers to Moments Applied by 

the Roll and Pitch Torque Generators 

In a manner parallel to that of Section 2. 3 the funda- 
mental equation of motion of a controller consisting of a gyro 
case, a gimbal, and a rotor is 




(Eq 6.2. 1) 



Since it is proposed to derive roll and pitch moments from the 
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control system by torquing the controllers about their vertical 
axis (y axis), it is convenient to sum torques in the coordinate 
frame of the case (gu frame). Accordingly, 



) M | =4r-| + W T ptt ^H| „ tt 

/ / casej q-jj dt _J qu J GU 

where H = H + H + H 
r g c 

From Appendix D, Equation D. 3. 5, 



W 



I, GU 



PCY + qSY 

-PST + qCY 
r + y 



The approximation made in Section 2. 6 is that 



XI | 

J GIM 

Now, in the gu frame 



0 

J j8 Ca 
P 

J p /3 So- 

Define J as the combined moment of inertia 

c gr 



H dGU = QGU ' GIM ” r J GIM 



0 

J jB 
p 

0 



(Eq 6.2.2) 



(Eq 6.2. 3) 



(Eq 6. 2. 4) 



(Eq 6. 2. 5) 



of the case. 
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gimbal and rotor. 



J = J + J + J 
cgr c g r 



(Eq 6. 2. 6) 



Combining the last two equations gives 



h Jgu - 



J ft Ca 
P 



J ftSa + J t 
p cgr 



(Eq 6. 2. 7) 



Differentiating with respect to time where it is assumed that 
angular momentum about the spin axis, Jp^ , is constant gives 






GU 



- Jpj3 Sn a 



4 J ft Ca a , T 
D f + J 



cgr 



(Eq 6. 2. 8) 



Then substituting the appropriate values into Eq 6. 2. 2 gives for 
the Number 1 Controller the following. 



V 

M 1 , 
„ v case 1 



- (R + y) J Cq'4-PSt ■ h qC7)(Jp(3So'l= J c 



g r 



GU 



-J ftSaa- (pCi + qS 7 )(J ft Sa* <T _ 7 ) 
P p ^'g 1 



+ J c gr ^ + + P 0 '!' + PS7) Jp^Ca 



(Eq 6. 2. 9) 
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The X and Y components of the torques in the last equation are 
normal to the single degree of freedom of the case, and there- 
fore are transmitted through the bearing to the vehicle. The Z 
component is aligned with the bearing axis, and the rotation of 
the case must be controlled by a torque generator about the 7 
axis. The single Z component can be written as 

M 7 easel = J c gr 7 + J p ^ ® + P C7 + qS7 Ca (Eq 6> 2 - 10) 

cage i consists of a torque generated part, a load, and a 
viscous friction. Applying these to equation 6. 2. 10 gives for 
each of the controllers the following. 

Controller 1 Let F, = J -77 + C 

1 egr dt 

M TG1 + M L1 = F 1 7 1 + J p 1 ;3 i + pC7 l + qS7 i )c »i 

(Eq 6. 2. 11) 



Controller 2 ( 7 ^ = 7^ + 180°) 

M TG2 + M L2 = F 2 7 2 + J p ^2 ^ a 2 " pC7 2 " qS7 2^ Co, 2 

(Eq 6. 2. 12) 

Controller 3 (7 3 = 7 - 90°) 

M TG3 + M L3 = F 3 7 3 + J p.A {Q 3 + pS7 3 ' qC7 3 )Cff 3 

J 

(Eq 6. 2. 13) 
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Controller 4 (y^ = 7 + 90°) 



M TG4 + M L4 " F 4 7 4 + J p 4 ^4 ^4 " pSt 4 + c l C7 4^ Ca 4 

(Eq 6. 2. 14) 



Combining Controllers 1 and 2 in accordance with the logic that 
7 = - 7 and o ' 1 = a , and making the further stipulations 

that 



M 



LI 



M 



TG2 



J 

p l 



h 




F 



1 



F 



2 



(Eq 6. 2. 15) 



which implies that the two controllers are geared together gives 



M 



7 1 = 



TGI 

^ 7 " 



2 J /SpCy.C 

E_ 1 



cv 



(Eq 6. 2. 16) 



• • 

The assumption that = a ^ is not exact because the a ^ and a ^ 

gimbal axes cannot be easily geared together, but if the variables 
are enclosed in an error position control loop the assumption is 
believed valid. Note that the gearing of controllers 1 and 2 com- 
pletely eliminates the cross coupling effects of q and a , and 
leaves only the primary attitude rate variable, p, of the loop. 
This is an important advantage of a system using two controllers 
operating back-to-back. 
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In arriving at the result of equation 6. 2. 16 is must 
not be overlooked that an important assumption is that the angles 
a i and a ^ are maintained approximately equal. This assumption 

gives the two- degree- of- freedom controllers the characteristics 
of a single- degree- of- freedom controller in that the a axis is 
approximately rigidly restrained unless all four of the controllers 
move in the same direction as a unit. In principle this may be 
accomplished by wiring the torque generators controlling the 
a axis in parallel and requiring the torque generators to have a 
high back emf. 

Combining controllers 3 and 4 under the conditions that 
the two are geared together and driven by controller 3 gives 
for similar assumptions to equation 6. 2. 15 that 



M 



TG3 + 2J eqCYgCa 



2F, 



(Eq 6. 2.17) 



Note that in this equation for pitch control that the cross coupling 
effects of roll and yaw do not appear in the equation. 

6. 3 Response of the Controllers to Moments Applied by 
the Yaw Torque Generators 

Following the procedures of the previous section 




(Eq 6. 3. 1) 



Then it may be written that 




Hi + W ^ H I 

GIM W I,GIM n J GIM 



(Eq 6. 3. 2) 



where H = H + H 
r g 
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Now from Appendix D the angular velocity is 



pCy + qSy + a 



W 



I, GIM 



(- pSy + qCyjCo- +• (r 4= 7 ) So- 



(Eq 6 . 3. 3) 



(pS 7 - qC 7 ) So- + (r + y)(Za 



Using the approximation of equation 6 . 2. 4 one would first 
guess that H I vanishes; however, such is not the case 

rg J GIM 

and we must use the following from Appendix E 



A 0 




0 



(Eq 6 . 3. 4) 



This step can be justified by considering the complete equations 
in Appendix E. 



Ym. , 
Lj g im 1 



U 



J r " - (pSy - qCy) Sa 4 = (r * 7) Co- J 
r g p 



0 



GIM 



£pCy + qSy + J j3 



J 

(Eq 6 . 3. 5) 



The Y component vanishes only for the assumptions made herein, 
which means there are no large torques transmitted about the 
spin axis of the rotor. Here we are interested only in the X 
component since the Z component will be transferred through 
the gimbal to the case where a component of it will appear at the 
7 axis. Note that the cosine a component can be identified in 
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the corresponding portion of equation 6. 2. 9. 



Unfortunately, there does not appear to be a simple 
method of mechanically connecting the gimbals of the four con- 
trollers so that each of these angles are exactly the same. 
Therefore, each of the gimbals must be handled separately. 
These may be written as follows 



j 

Controller 1 Let = J r g ^ + C 

. m tgi * j p^ (pS7 i -qcysfri + + 



(Eq 6 . 3. 6 ) 



Controller 2 



= 



M TG2 + J p 13 (pSTj + qC7 1 )Sa 2 + (r - 7^0*2 

L — 



(Eq 6. 3. 7) 



Controller 3 



' M TG3 + J P^ ^ pCT 3 ~ + (r + ^3) 

0,3 ” G s 

(Eq 6. 3. 8) 

Controller 4 

* 

* M TG4 f J P /3 (_ P Ct 3 ' qS 7 3 ) Sn 4 + (r " ^ 3 ) Ca, 

a 4 GJ 

(Eq 6. 3. 9) 
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In the above four equation it is assumed that 7 g = - 7 -^ ’ ^4 = ~^3 

and that J 8 is held constant. 

P 

If it is assumed momentarily that the a gimbal axes 
can be mechanically geared together then one may sum the yaw 
moments of the four controllers. Since the angular momentum 
of each controller is identical a summation of yaw moments can 
be accomplished by first finding the average gimbal rate by sum- 
ming the four gimbal rates. This gives the following equation. 

* 2M TG * J p^C (pST l " qS T 3 > Sa + 2r C a] 
a (Eq 6 - 3 - 10 ' 

In this expression the moment applied by the torque generator, 

, is applied to each of the four gimbal axes. Note that 

« » 

unlike the previous equations for 7 ^ and the cross- coupling 

rates p and q do not vanish although they are multiplied by the 
product of two angles which may remain small. However, the 
large cosine components of the cross- coupling rates do drop 
from the equation as do the terms in 7 ^ and y^ . 

It is proposed to gear the a gimbals together electrically 
by finding the average gimbal angle and comparing this with the 
angle of each gimbal separately in an error feedback loop. The 
equations for this operation are contained in the next section. 
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Torque Generator 



6. 4 

The control power requirements for the gyro controllers 
are insignificant compared to the power required to drive the 
rotors at operating speed as will be shown, and since the peak 
torque requirements are small, either a direct current motor 
or a two phase alternating current motor is applicable. The 
ideal torque-speed characteristics of both types are the same 
to a first approximation. 

J M 0 + B9 + M l = K t V m (Eq. 6.4. 1) 

where 9 = motor position angle 

Jjyj = motor moment of inertia 

B = viscous friction constant 

Mj^ = load torque 

Krp = torque constant 

Vjyj = input control voltage 

The motor is geared to the controller with a gear 
ratio p > 1; therefore, the torque applied to the controller 
gimbal is 

M TG = pM L (Eq. 6.4.2) 

The angular relation between the motor position angle 
and gimbal position angle is 

9 = pa for yaw control 
9 = pY for roll and pitch control 



This gives for yaw control 



(Eq. 6. 4. 3) 



and for roll and pitch control 



M TG = K Te V M-( J M'’ 2 “ +B ‘ >2 - ? ) 



(Eq. 6. 4. 5) 



Combining equations 6. 2. 16 and 6. 4. 5 gives for roll control 

(Eq. 6. 4. 6) 



T l = 



K T p V M ’ 2 J p 0 P Cy i Ca 



It 



where 



F lY< 2J Cgr +J MP 2 )jT + (2C + B p 2 ) 



Likewise for pitch control 

, K T p V M + 2Jp^qC7 3 Ccv 

t 3 = FT 

3t 



(Eq. 6. 4. 7) 



For yaw control the four controllers give separately 
the following relations. 

Let G l.=( J rg +J MP 2 >3T + ( C+B 0 2 ) 

• k t’ v m +j p^[(p s 7^ c 7>i + ( r+ ^i) Cff i] 
a i tr 

la 

(Eq. 6. 4. 8) 

, K T p V m + J p |3 [(pSTj + qCy^S^ + (r-y^Co^] 

* 2 = s; 

(Eq. 6. 4. 9) 
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K T P V M + J P^ LC pC7 3 " qST 3^ Sc>, 3 + ( r +7 3 ^) Cq, 3] 

G 3q' 

(Eq. 6. 4. 10) 

K T pV M +J pH(' pCT 3' qS7 3K + ( r '^ Ca 4] 

G 4q' 

(Eq. 6. 4. 11) 

In the above equations subscripts have been omitted, but each of 
the equations requires its own characteristic parameters. 

Since all controllers are slaved to the position angle of 
the average controller, the relation for the motor voltage of a 
typical controller is as follows. 

V M = V MZ- Kfl ‘ , (‘V“ave) (Eq- <S. 4 . 12 ) 
where = error voltage and K0 = feedback constant. 

Substituting equation 6. 4. 12 into the above four equations gives 
the following when written in terms of position angles. 

Controller 1 

a _ K T pV MZ + K 9 p ^ave + J P^K ST l" qCT l!) SQ 'l + (5’ atl,Y l^) CQ 'l J 

1 tfCO + K e p 

(Eq. 6. 4. 13) 




Controller 2 



ot n - 



K T^ V MZ tK 9' Ia ave' fJ P^(P S ^ + q C n>‘ 1 '2< r -ii> g 2] 

4t ( G 2«) * K 9 P 



(Eq. 6. 4. 14) 
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Controller 3 



Q'o = 



K tP V MZ+ K QP‘W ^ J pg[(p c T 3 -q s T :i >»3<rH.T , 3>a'3] 



^rC G 3a) + K 0 p 



(Eq. 6. 4. 15) 



Controller 4 



O'. - 



K T pV MZ* K 9 p,> a ve' |,J P^K: pC ' l '3' <lS ' l '3) S,> 4' t, ( r ~ , >3) C:a '4] 



7 



d_ r 
dt u 4» 



+ K fl P 



(Eq. 6. 4. 16) 



If we assume that the gimbal angles are held approximately 
equal such that Ca^ = Ca^ = Ca^ = Ca ^ , then the operation of 
finding the control moments contributed by the four controllers 
can be accomplished by first finding the average value of the 
gimbal rates by combining the last four equations. This gives 
for the average gimbal rate the following 



2K T pV MZ +Jp^[(pS7 1 -qS7 3 }s^ + 2rCa'] 



2 G 



a 



(Eq. 6. 4. 17) 



Where G 2 = (j rg + J M P 2 ) £ + (c + B p 2 ) 
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Attitude Sensors 



6. 5 

The primary source of guidance, navigation, and attitude 
information for a spacecraft is from external sightings. As was 
shown in Chapter 4 the ambient fields around a spacecraft are 
too feable to drive a practical sensor, and this implies that some 
type of star tracker must be provided if the spacecraft is to be 
self sufficient. A control system can control the attitude of a 
spacecraft only to the precision provided by the attitude sensor, 
and the design of a highly linear sensor is a problem within it- 
self. As done by other investigations of hardware for applying 
torques to a spacecraft, this thesis assumes a linear sensor 
for each of the three attitude reference axes of the spacecraft. 
Accordingly the following equation is assumed to represent 
the sensor- 



K K 
sx ex 


0 


0 




</> r - </> 




V 

cx 


0 


K K 
sy ey 


0 




o - e 

r 


= 


V 

cy 


0 


0 


K K 
sz ez 








V 

cz 


— 




— 




__ 







(Eq 6. 5. 1) 



K g = sensitivity of sensor 

= sensitivity of error amplifier 
<j> r = reference angle 
$ = vehicle attitude angle 

V c = voltage signal to tandem compensation 

Any practical sensor has linear characteristics over 
only a small range, and the effects of saturation will be consid- 
ered in a later section concerning the position control mode. 
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6. 6 Tandem Compensation 

To accommodate the provision for shaping the error 
signal before it is fed to the control system torque generators, 
a transfer function which is presently undefined is inserted into 
each control loop. The following equation applies. 



" c x 


0 


0 




~ v cx 




v 

MX 


0 


C Y 


0 




v cz 


= 


V MY 


0 


0 


c z_ 




v cz 




v 

MZ 



(Eq. 6. 6. 1) 

V c = voltage signal to tandem compensation 
C = transfer function of compensation 
= voltage signal to torque generators 

It is not anticipated that any tandem compensation will 
be required in the examples of this thesis. In an actual problem 
certain regions of the s-plane may be denied to the control sys- 
tem designer because of unstable modes of the airframe, etc. , 
and for this reason space in the equations has been provided 
for tandem compensation. 
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Position Control Equations 
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For convenience we write equation 6.8. 1 in shorthand notation 
as follows. 



<t> r - 4> 
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e - e 
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+ 4 H 2 [b] 
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4^ - ^ 




r 



2 / M (l2-34-l234)J A 2 H 1X) 



(Eq 6. 8.2) 

where ^ A J and are the indicated matrices and (a^) and 

(b^) are typical elements of these matrices. 

To obtain the equation for the complete closed loop 
system we substitute equations 6. 8. 2 and 2. 4. 5 into equation 
2. 3. 3 which gives the following. 
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(Eq 6. 8. 3) 
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Equation 6.8.2 represents the complete loop of the 



attitude control system, and as can be seen there is non-linear 
interaxis coupling between the three loops as well as non-linear 
terms arising from the Euler inertia-cross coupling. If it is 
assumed that compensation is applied to the loops as set forth 
in section' 3. 6, and also compensation for the inertia-cross 
coupling terms, the following expression is found for the roll 
attitude angle. 



Similar expressions can be written for 6 and i|j . 

If the angular motions of the a axes of the controllers 
are restricted to small angles then the following equations apply. 




dt 



I F. 
X 




(Eq 6. 8.5) 




(Eq 6. 8.6) 



I G ^ + (2 HCn)^4£ + 4 HQ* K i)j = GaM +4 HCa K ib 
z a dt z z^r 



dt 



(Eq 6. 8. 7) 
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where 




F 



37 



G 

a 



< 2 J cgr +J M' >2> ^ 



< 2J cgr +J M 



n + t 2*_d_ + 
J gr J M p ) dt 



+ (2C + Bp 2 ) 

+ (2C + Bp 2 ) 
(C + Bp 2 ) 



the constants for the different axes may not be equal. 

Block diagram representation for the above equations 
are contained in Figure 6. 8. 1 and Figure 6. 8. 2. Figure 6. 8. 1 
contains the roll and pitch channels whereas the yaw channel 
is shown in Figure 6.8.2 in two parts. Part (a) shows the four 
individual controllers driven by the same error voltage. To 
keep these controllers to the same angle their output angles are 
averaged and then compared with the actual output gimbal angle. 
The error is then fed back to the gimbal torque motor. Part 
(b) of Figure 6.8.2 is the result of ideal controllers and repre- 
sents equation 6. 8. 7. 
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Roll Equation and the Pitch Equation for Uncoupled Motion. 
Equations 6.8.5 and 6.8.6. 
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PART (o> INDIVIDUAL CONTROLLERS AND AVERAGING CIRCUIT SHOWN 



^ext 




PART (b) SIMPLIFIED REPRESENTATION OF EQUATION 6.8.7 



Figure 6.8.2 Block Diagram Representation of the Position 
Control Yaw Equation for Uncoupled Motion. 
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Rate Control Equations 



&JL 

The purpose of the rate control is to provide the 
astronaut with a means of slewing the vehicle at a particular 
attitude rate. Since the position of the gimbal angles are 
directly related to the amount of control system angular mom- 
entum transferred to the vehicle, the gimbal angles represent 
attitude rates. Therefore the rate control is essentially an 
error feedback control loop to control the gimbal angle, and 
the block diagram of the roll channel is shown in Figure 6. 9. 1. 
During this operation the error signals from the vehicle attitude 
sensors are open circuited. 

During the rate control maneuver the voltages to the 
torque motors are given by the following equation. 
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(Eq 6. 9. 1) 



The voltages given by equation 6.9.1 may be substituted into 
equations 6. 4. 6, 6. 4. 7, and 6. 4. 17 to find an expression for 
the gimbal angles as a function of the rate variables and the 
rate inputs. For shortness we write 



K = K K p 
1 yl tx'x 



K 9 = K Q K p 
2 y3 ty K y 



K„ = K K, p 
3 a ta a 



which gives the following matrix equation. 
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(Eq 6. 9. 2) 
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The above matrix equation must be solved for the gimbal angle 
rates so that these rates can be applied to the control system 
represented by equation G.2. 11. Because of the non-linearities 
caused by the sine and cosine terms of equation 6. 9. 2, the 
Laplace techniques are not simple to apply unless the angles can 
be considered constant. Therefore, in the equations that follow 
the notation employs the letter s so that a closed form solution 
can be written, but this symbol s refers to the Laplacian oper- 
ator only in the special case where sine and cosine values of 
the gimbal angles can be considered constant. Otherwise, for 
cases where the gimbal angles make large changes, the symbol 
s must be interpreted as a differentiation with respect to time. 

Solving the matrix equation 6.9.2 for the gimbal angle 
rates and substituting in equation G. 2. 11 gives the following 
expression for the moments generated by the spacecraft control 
system. 
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As in the previous solution, rewrite equation 6. 9. 3 as follows. 
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(Eq 6. 



Substituting equation 6. 9. 3 and 2. 4. 5 into equation 2. 3. 3 
the following equation for the response to rate control. 
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The uncoupled set of equations for roll, pitch and yaw 
is given by the following. 



M (SF. +K.) + 7 . (2HCt 1 K 1 S) 
x I7 1 ' lr ' 1 1 

S (I F. S + I K + (2 HCyJ 2 ) 
x I 7 x 1 1 / 



(Eq 6 . 9. 6 ) 



M y <SF 3t + K 2> - V 3 r (2 ”^3 Kg S _> 
S(I y F 3Y S + + < 2HC V 2 ) 



(Eq 6 . 9. 7) 



M (SG + K„) - a (4 HCa K Q S) 
z a 3 3 

SAG S + I L + (2HCyJ 2 ) 

Z ct Z o Of 



(Eq 6 . 9. 8 ) 



The following equations are valid for the rate input 
commands which drive the gimbal angles to large values in 
those cases where there is negligible inter-axial coupling. 



For Roll 
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(Eq 6 . 9. 9) 



For Pitch 
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(Eq 6 . 9. 10) 



For Yaw 
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(Eq 6 . 9. 11) 
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For the completely general case it is necessary to 
solve equations 6. 9. 2, G.2.11, 2.4.5, and 2. 3. 3 by machine 
computation. 




Figure 6.9.1 Block Diagram Representation for Roll Channel 
of the Uncoupled Rate Control Equation. 
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6, 10 Summary 

This chapter has served to derive the closed loop 
equations of motion of the spacecraft. Examination of the 
closed loop position control equation 6. 8. 3 and the closed 
loop rate control equation 6. 9. 5 confirms the fact that the 
equations are coupled from control interaction as well as the 
Euler rate-product terms. Therefore, a completely general 
solution requires assistance from a computer. One of the 
reasons for choosing the four controller system is that many 
of the large cross- coupling terms in the primary control 
matrix vanish, and since the four controller system is a zero 
momentum type system the cross- coupling terms of the gyro- 
scopic coupling matrix are small. Cross- coupling is a problem 
in any type of flying machine if we expect to actuate the con- 
trol of two axes to their full travel simultaneously. Generally 
speaking vehicles are not operated in this manner, and if 
crossed controls are required one of the axes will predominate. 

In the following chapter the equations of motion are 
solved to determine the response of the vehicle. This is ac- 
complished by assuming no interaxial coupling, and solving 
the equations analytically. Making this assumption, many of 
the equations can be solved in closed form whereas others 
require graphical methods. It is believed that this approach 
gives a better understanding and presentation of the system 
than a pure computer solution which gives graphical response 
of the vehicle to input disturbances. 
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CHAPTER 7 



DYNAMICS OF THE SPACECRAFT 
7. 1 Introduction 

This chapter evaluates the response of the spacecraft 
using the equations developed in the preceding pages. This is ac- 
complished by assuming that the interaxial coupling among the 
three control axes can be neglected. This assumption is valid 
provided the control system is not appreciably saturated, and 
that the astronauts do not attempt to simultaneously actuate the 
controls to yield large rates about two axes. In any event, a 
necessary condition is that the attitude control system must per- 
form satisfactorily about each of the individual control axes so 
that initial analysis of the equations of motion on an uncoupled 
basis is considered reasonable. Also, several parameters of the 
system must be numerically chosen, and their selection is greatly 
simplified using uncoupled equations. These parameters include 
the damping coefficient of the controllers, the gain of the rate 
control loop, and the gain of the position control loop. It is visu- 
alized that the spacecraft will be operated in one of three modes. 

1. Zero Input Mode 

2. Rate Control Mode 

3. Position Control Mode 

In addition to these three modes, a mode called the Adaptive 
Mode will be discussed. 

The zero input mode is that which results when no signal 
is provided to the torque motors of the controllers. The control 
system then operates as a rate stabilizer, and the gyro controllers 
provide inherent sensing. The stability of the zero input mode 
will serve as a criterion for the selection of the damping coefficient 
of the controllers. The damping coefficient is then held constant 
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for the other two modes. The purpose of the rate control is to 
command the spacecraft to perform a particular attitude rate. 

The stability of the rate control mode will be used to ascertain 
the gain of the rate control loop, and the stability of the position 
control mode will determine the gain of this loop. The adaptive 
mode is the vehicle response to a controller failure. 

The damping chosen for the controllers differs appreciably 
from that used in systems which overdamp the gyro to the point 
that it can be approximated as a first order system. In the present 
design the controllers act as a true second order system with 
two poles having both real and imaginary parts. This of course 
complicates the solution of the equations, some of which cannot 
be solved in closed form. Therefore, in many of the equations 
which follow it is necessary to approximate the second order 
system by a first order system to arrive at useful results. 

As was done in Chapter 5, the analysis will be accom- 
plished using the equations for roll. The equations for pitch con- 
trol are almost precisely the same as those for roll, and generally 
the same parameters are required except that the numbers will 
differ because of different moments of inertia of the spacecraft 
about the roll and pitch axes. The yaw equations differ slightly 
from those of roll, and a summary of the differences will be 
made. 
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7. 2 



Zero Input Mode 



The uncoupled equations for the zero input motion are 
given by equations 6. 8. 5, 6. 8. 6, and 6. 8. 7 where the position 
input terms have been dropped, and these equations are written 
as follows. 
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(Eq. 7. 2. 1) 
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(Eq. 7. 2. 2) 
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(Eq. 7. 2. 3) 



The stability of the spacecraft for the complete range of 
gimbal angles can be investigated by assuming small perturbances 
wherein the gimbal angle can be taken as a constant. Consider the 
equation for roll motion and let which is actually equal to 

= (2 J c g r * Jjyj p 2 ^) s + (2 C + B p 2 ^) be written, for short- 
ness, simply as F^ = J s + k. This gives 



P =■ 
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(Eq. 7. 2. 4) 
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The motion follows that of a second order system with a damping 
ratio 



C = 



k 

4HCy 1 




(Eq. 7. 2. 5). 



and a natural frequency. 
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(Eq. 7. 2. 5) 



The effect of the disturbance is least when the gimbal angle is 
small for which the following are defined 
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(Eq. 7. 2. 6) 
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(Eq. 7. 2. 7) 



The solution for the response of the system to an impulsive mo- 
ment, M q) for the lightly damped controller is given by the fol- 
lowing equation. 
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(Eq. 7. 2. 8) 



Examination of this equation shows that the roll rate is always 
stable for any gimbal angle less than 90 degrees because the ex- 
ponential damping factor is not a function of the gimbal angle. 
Critical damping occurs when Cy^ = and the solution for this 
motion is given by the following equation. 
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(Eq. 7. 2. 8) 
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The highly damped motion can be approximated by assuming that 
the gimbal moment of inertia is negligibly small giving the fol- 
lowing equation for roll rate. 

M 1 

4H 2 (CYi) 2 ! 

J x k J 

(Eq. 7. 2. 9) 

The solution to this equation for an impulsive moment disturbance 
is given by 

4 H 2 Ct 2 

M - _ ^ T 1 t 
p = -y Q " e Ik 

I x 

x (Eq. 7. 2. 10) 

In terms of previously defined parameters this motion can be 
written as follows. 

- u> 

S- t 

M 2 ^ 

p = -j- 2 - e (Eq. 7. 2. 11) 

x 

The analysis of the response due to a step moment dis- 
turbance is much more difficult to analyze because the gimbal 
angle cannot be assumed constant. As a beginning, however, let 
us assume that the moment is small enough that first order dyna- 
mics can be used in which case equation 7. 2. 9 is considered. 
Further, let the gimbal angle remain small so that the cosine of 
the gimbal angle can be approximated as a value of one. A step 
of at time zero gives the following. 
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(Eq. 7. 2. 12) 
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The final value is seen to be independent of the moment of inertia 
of the vehicle. The roll rate of a vehicle without an attitude control 
system in response to a step disturbance moment is given by 



P = 




(Eq. 7. 2. 13) 



The roll rate of a spacecraft without a control system compared 
to that with gyro controllers both of which experience a step mo- 
ment disturbance is given by dividing equation 7. 2. 13 by 7. 2. 12 
as follows. 



^without control system _ M o 
^with gyro controllers ^ 0 




w 



The factor 



2 t 



(Eq. 7. 2. 14) 

is the reciprocal of the time constant for the 



gyro controller and the factor is large for fast controllers, so 
that the ratio is plotted initially as a steep curve in Figure 7. 2.1. 
The effect of saturation of the gyro controllers is sketched in to 
show that the ratio approaches unity as time increases. 
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SYSTEM 
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Figure 7.2.1 A plot of the ratio of roll rate of an uncontrolled 
spacecraft compared to a gyro controlled spacecraft both of 
which experience a step moment disturbance. Dotted curve is 
undocumented. 
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The solution for a step moment disturbance wherein the 
controller gimbal angles move through their complete range is 
solved by a phase plane analysis, however, first consider the 
special case of recovery from a large roll rate in the absence 
of external moments. The solution for the gimbal angle has been 
determined as follows. 



y^ - arc tan 




tan y 

o > 



(Eq. 7. 2.15) 



Using this relation for the gimbal angle gives the following equa- 
tion for the roll rate. 
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These equations are plotted in Figure 7. 2. 3 for an initial 
angle of 89 degrees. 



7. 2. 16) 
gimbal 




Figure 7. 2. 3 Plot of gimbal angle, y^, and non-dimensional 
roll rate, p/p , for spacecraft controlled by two gyro control- 
lers with initial condition of y^(0) = 89 degrees. Total angular- 
momentum of system is zero and there are no external moments. 
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This special case of recovery from a large roll rate in 
the absence of external moments also has a simple phase plane 
solution as follows. As shown in Figure 7. 2. 4, plot the following 
equations for selected values of gimbal angle. 

. (2HCt,) 2 

p * — ^ — p = 0 (Eq. 7.2.17) 

x 

if," p = p Sy, (Eq. 7. 2. 18) 

The trajectory is determined by intersections of lines of the same 
gimbal angle. The times of the trajectory can be found from the 
ratio Ap/p avera g e , and from this the curve roll rate versus 
time can be plotted. Thereafter, roll rate can be integrated to 
find the roll displacement as follows. 

C * 

</> = \ pdt (Eq. 7. 2. 19) 
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Figure 7.2.5 The relationship between gimbal angle and roll 
displacement for the zero input mode with zero initial conditions. 
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A closed form solution of equation 7. 2 . 1 is difficult even 
when only the first order dynamics are considered. There is, 
however, a direct relationship between the vehicle attitude angu- 
lar displacement and the gimbal angle. For example, the follow- 
ing equation is always true for zero initial conditions and for zero 
input to the controller. 

Jy* + ky= -2HCy p (Eq. 7. 2. 20) 



Considering only first order dynamics the equation can be inte- 
grated to give the following result. 



7=2 arc tan* 



2 H 



'-y (Eq. 7. 2.21) 



This equation is plotted in Figure 7. 2 . 5. 



An approximate phase plane analysis can be made of the 
position input mode operating to minimize a step torque distur- 
bance by assuming that y is constant. The first order differential 
equation approximation of equation 7. 2. 1 is given by 



, 4H 2 (C 7l ) 2 M 

p+ l ,. k 1 P =Y JL (Eq. 7.2.22) 



This equation is plotted in Figure 7. 2. 6 for various contours of 
7 j which also represent contours of constant time since 

A 7i 

At =-r~ L - (Eq. 7. 2. 23) 

7 1 



Since increased gimbal angles demand increased values of the 
roll rate, p, it is clear that p must also increase. At the end of 
the transient when p = M x k/4H 2 , p is very small, but its 
average value over the interval from 7 ^ = 0 ° to 7 ^ = - 10 ° can 
be found by reading from Figure 7. 2 . 6 the value of 

- AE 

M x k 
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from the relation 



and computing p 



ave 



Pave 'If (Eq. 7.2.24) 

Figure 1 . 2.1 is an analog simulation of the zero input mode and 
shows that 7 holds constant at -M/2H until the gimbal angle gets 
beyond about 40 degrees. A phase trajectory determined by the 
method outlined above is shown in Figure 7. 2. 6. 
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Figure 7.2.6 Estimated Phase Trajectory for Zero Input 
Mode Response to a Step Moment Disturbance. 
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All of the equations for zero input control have been presented 
in parametric form with no actual numerical values substituted in 
the equations. Therefore, the equations are valid for any general 
twin gyro controller system where interaxial coupling is negligible. 
The zero input mode is purely a rate stabilization system. It has 
been shown that an impulsive torque disturbance acting on the 
vehicle will cause the spacecraft to roll through a displacement 
angle, but that the roll rate will damp to zero. A step moment 
disturbance will cause a steady state roll and consequently an in- 
creasing roll displacement. The zero input system is primarily 
a rate stabilization system and the parameter which can conven- 
iently be chosen to give the desired response is the damping 
coefficient, k, of the gyro controller. It is not necessary at this 
point to choose a particular value of k other than to say that this 
quantity is chosen on the basis of equation 7. 2. 6 which solved for 
the damping coefficient is written as follows. 



From a practical standpoint all of the quantities in equation 
7. 2. 22 except are dictated in the design by other considerations. 
I concerns the overall vehicle design and is determined by the 
mission and the detail design of the spacecraft. The angular mo- 
mentum of the controller, H, is determined from controllability 
requirements, and the combined inertia of the gimbal, case, and 
rotor, J, is to be determined by optimum design of the controller 
to minimize the ratio of the total mass of the controller compared 
with the angular momentum of the controller. This leaves the 
parameter which is to be chosen by the desired location of the 
closed loop poles of the zero input mode. Although the damping 
coefficient, k, is chosen primarily on the basis of the desired 
response to disturbances of the zero input mode, the choice does 
affect the other modes. For example, the closed loop poles of 
the zero input mode will be the open loop poles of the position 




x 



(Eq. 7. 2. 27) 
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control mode. Therefore, let us delay the numerical choice of 
the damping ratio, C, , until the equations of the other modes 
have been examined. 







Figure 7.2.7 Spacecraft Response in Roll to a Step 
Moment Disturbance of 10 lbrft. I = 10^ lb-ft-sec^, 

J = 10 lb-ft-sec^ H = 10 4 lb.-ft-sec, £ Q = 0.866. Data from 
Analog Simulation. 
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7. 3 Rate Control Mode 



The equations for rate input control given by Equation 6. 9. 9 
can be solved graphically by a plot similar to 7. 2. 4 where the 
trajectory is defined by the intersections of lines of equal gimbal 
angle. For this solution equation 6. 9. 9 is written in a slightly 
different form which, neglecting the external moment disturbance, 
is as follows. 
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(2HCt)‘ 
I k 



p= ( T lr' 1 'l) 



K 



1 2 H Cy 
I k 



(Eq. 7. 3. 1) 



This equation is plotted for p versus p in Figure 7. 3. 1 and gives 
contours of constant gimbal angle which slope from left to right. 
The second equation needed is that of equation 7. 2. 18 which is 
the requirement that angular momentum is conserved. This 
equation is plotted in Figure 7. 3. 1 as vertical contours of con- 
stant gimbal angle. The intersections provide the solution. The 
gain of the gimbal position control loop, , has been allowed 
to have a relation to other parameters of the vehicle as follows. 

2 

_ m ±jl (Eq. 7 3 . 2 ) 

• * x 

The factor m may be called a gain ratio, and Figure 7. 3. 1 indi- 
cates that for high gain ratios the roll rate is an exponential 
represented by a constant slope on the plot of p versus p. The 
effect of the gain ratio factor, m, is shown more clearly in 
Figure 7. 3. 2 which is a plot of the steady state gimbal angle 
versus the value of the step input. This plot shows that for very 
high gain the gimbal angle is linear with the input angle, but for 
a gain ratio of one the variation is highly non-linear. If the gain 
ratio is less than one, the position control system does not have 
the capability of holding the complete range of gimbal angles in 
a stable manner. Figures 7. 3. 3 and 7. 3. 4 further illustrate the 
effect of low gain ratios. Figure 7. 3. 3 shows a plot of steady 
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Figure 7.3.1 Phase trajectory for spacecraft response 
to a full step displacement rate control input. The external 
moment disturbance is zero. The response is defined by 
the intersections of the lines of constant gimbal angles 
plotted for the two equations of motion shown on the 
figure . 
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Figure 7.3.2 Plot of steady state gimbal angle, 7 , , 
resulting from step input, 7 ^ , for rate input mode. 
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Figure 7.3.3 Steady state values of roll gimbal angle for 
gain ratios less than one in the gimbal angle position control 
loop used in rate input control. 
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Figure 7.3.4 Plot of roll capability for gain ratios less 
than one in the gimbal angle position control loop used in 
rate input control mode. 
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state gimbal angle for a step input of ninety degrees versus the full 
range of gain ratios less than one. Careful inspection of the cur- 
vature shows an inflection point at 45 degrees which was shown 
to be a peak control power point in section 5. 31. Figure 7. 3. 4 
shows the roll capability for a step input of ninety degrees for 
the range of gain ratios less than one. 

From a stability standpoint the rate input mode has a natural 
frequency from equation 6. 9. 6 equal to 



This natural frequency is greater than that for zero input control 
given by equation 7. 2. 5. 

The product of the damping ratio and the natural frequency 
are the same for zero input and rate input. Thus 



If the damping coefficient, k, is chosen on the basis of zero input 
control then according to equation 7. 2. 6, 




(Eq. 7. 3. 3) 




(Eq. 7. 3. 4) 




(Eq. 7. 3. 5) 



x 



Substituting this k into equation 7. 3. 4 gives 



4H£, 



o 




(Eq. 7. 3. 6) 



Then further substitution of in equation 7. 3. 2 into equation 
7. 3. 3 gives 



Therefore, 



C 




(Eq. 7. 3. 8) 



which is to say that the damping ratio for rate control is propor- 
tional to the damping ratio for zero input control by the factor 



Actually, large values of the gain factor m gives more precise 
positioning of the gimbal angle as shown by Figure 7. 3. 2, but the 
greater this gain the more lightly damped (more oscillatory) be- 
comes the rate control mode. It appears that an acceptable solu- 
tion lies in the choice of the gain ratio, m, as unity. For a gain 
ratio of unity the gimbal angle can be positioned at the full ninety 
degree position, yet for the range of 45 degrees or less the curve 
for m = 1 of Figure 7. 3. 2 is approximately linear. A gimbal 
angle of 45 degrees gives more than 70 of the maximum roll rate; 
therefore, it is not likely that the vehicle will be operated at 
angles of more than 45 degrees except for the full roll rate posi- 
tion of 90 degrees. For unity gain ratio 



By cross -plotting the steady state response to an impulsive 
disturbance torque on the spacecraft shown in Figure 7.3.5 the 
roll rates can be compared with those for zero input mode and 
those for no control system aboard the vehicle. For unity gain 
the rates are approximately one-half that of a vehicle without a 
control system. 



1 H m + (Cy)^ . The greatest difference between the two damping 
ratios occurs at zero gimbal angle. Hence, 




(Eq. 7. 3. 9) 




(Eq. 7. 3. 10) 



195 




Figure 7.3.5 Plot of steady state roll rate versus 
magnitude of torque disturbance impulse for rate control 
mode. Infinite gain ratio corresponds to controllers that 
are held rigidly fixed in the vehicle, and is the same as 
vehicle response with out a control system. Zero gain 
ratio corresponds to zero input mode of section 7.2. 
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7. 4 Position Control Mode 

The position control mode is formed by closing an attitude 
angle feedback loop around the zero input mode. A convenient 
method for studying the position control mode is a locus of roots 
which is a useful technique found in most control theory texts. 
Using root locus theory it is known that the closed loop poles of 
the zero input mode form the open loop poles of the position con- 
trol mode. Therefore, the open loop equation for the position 
control mode can be written for roll motion as 



( 0 ) 



open 

loop 



(k x 2HC Ti ) i r+ F 1t M x 

s(i x f i t s + 4 h 2 (cti) 2 ) 



(Eq. 7. 4. 1) 



where the gimbal angle must be considered approximately con- 
stant. The closed loop equation is given by 

( K x 2H CT 0 tr + F 1t M x 
" 0 x F i T s2 + 4H 2 (CTi) 2s + K x 2 H C 7i) 

(Eq. 7. 4. 2) 

. ' * 

The poles of the denominator of equation 7. 4. 1 are to be deter- 
mined so that the position control mode has good response to input 
commands and to minimize the disturbances of external torques 
on the spacecraft. An interesting theorem which applies to a con- 
trol system with no zeros and three open loop poles, one of which 
is at the origin, is as follows: the natural frequency of a unity 

feedback control system having no zeros and three open loop 
poles, one of which is at the origin, can never be greater than the 
natural frequency of the system operating open loop. This theorem 
is easily proven by substituting jco in the characteristic equation 
of the closed loop system and equating the imaginary part to zero. 
The theorem is illustrated in Figure 7. 4. 1 by noting that all stable 
poles of the closed loop system lie within the circle of radius oo 
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Figure 7.4,1 Illustration of four possible geometric patterns 
for the locus of roots for the position control mode. The three 
poles shown are the poles of the closed loop zero input mode 
which constitute the open loop poles of the position control 
mode. The closed loop poles of the position control mode will 
lie on the root loci depending on the value of the gain of the 
position control loop. 
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which is the natural frequency of the open loop system. The pri- 
mary purpose in Figure 7. 4. 1 is to illustrate the four possible 
geometric patterns that the locus of roots of the position control 
system can have. If the damping ratio of the zero input mode, Z, , 
is greater than one, then Figure 7. 4. 1 (a) applies. If the damping 
ratio is less than one but greater than 0. 866 then Figure 7. 4. 1 (b) 
applies. If less than 0. 866 then Figure 7. 4. 1 (c) represents the 
locus of roots. The case where t, is exactly equal to 0. 866 is 
shown in Figure 7. 4. 1 (d). The overdamped case (a) has poor 
response at low gains, and it becomes too oscillatory at the high- 
er gains. Values of t, less than 0. 7 are not satisfactory because 
the rate control mode and the position control mode both become 
too oscillatory. Values of £ between 0. 7 and 1 . 0 give an accept- 
able solution so that the vehicle has good response to input com- 
mands as well as good response to damp torque disturbances. 

As an example of this section let us solve the equations of motion 
using a damping ratio of C, = 0. 866 such that the locus of roots 
of Figure 7. 4. 1 (d) applies. This is not suggested as an optimum 
choice but is considered to be an example of a satisfactory loca- 
tion of the closed loop poles. Choosing a value of K such that all 

X 

three closed loop poles coalesce gives the following equations. 



K 

x 



0. 77 



_ni 

Vi j 

x 



(Eq. 7. 4. 3) 



0. 192 u> 3 CYi <i> + (S + 3(6. 577 oj ") ") M /I 
^ _ o 1 r V y 0/ y xx 

S 3 + 3(0. 577 w^S 2 + 3(6. 577w Q ) 2 (C7 1 ) 2 S+0. 192 w 3 CTj 

(Eq. 7. 4. 4) 

where is natural frequency of zero input mode given 
by equation 7. 2. 7. 
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To further explain equation 7. 4. 4, if the cubic denominator 
is written in the following form 

(S 2 + 2t pUn S+ . n 2 ) (S+ §) = 0 (Eq. 7. 4. 5) 

then the solution chosen as an example has the following values 
for the parameters. 



I = 1 (Eq. 7. 4. 6) 

P 

gj = 0. 577 gj (Eq. 7. 4. 7) 

no M 

5 = 0. 577 co (Eq. 7. 4. 8) 

o 

The primary purpose of the position control mode is to enable the 
spacecraft to track a reference line. This reference line may be 
the line of sight to a star, and usually in this mode the vehicle 
attitude rates are either small or are constant. It is considered 
that equation 7. 4. 4 can be evaluated for constant gimbal angles 
if we assume that the error signal from the sensor is limited to 
some maximum value. In actual practice this limiting will be 
accomplished by the saturation of the sensor. If the maximum 
error signal from the sensor is <j> then the steady state roll rate 
resulting from the constant sensor error can be derived as the 
following equation. 

v-Wf 

x 

lEq. 7. 4. 9) 

This curve is shown in Figure 7. 4. 2 and shows that for a given 
error signal, <j> , there are two solutions for roll rate. The 
solution at the higher rate is unstable for the same reasons that 
the initial point at p is unstable in Figure 7. 2. 4. Therefore, 
to avoid all possibilities of the position control mode being un- 
stable the error signal from the sensor is limited to that giving 
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a maximum gimbal angle of 30° which allows the position control 
loop to drive the rate of the spacecraft up to one-half P max - This 
limiting also enables the stability analysis of the position control 
equations to be accomplished with a constant gimbal angle. Thus 
for small gimbal angles the position control equation 7. 4. 4 can 
be written as follows. 



0. 192 W 3 0 + 
o r 


— — 
S +3(6. 577 Wa ) 


M /I 

X 


( S + 0. 577 





(Eq. 7. 4. 10) 



For zero sensor input the response to an impulsive moment dis- 
turbance is given by 



0 = 
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(Eq. 7. 4. 11) 

The response to a step moment disturbance is given by 
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(Eq. 7. 4. 12) 



For zero moment disturbance the response to a step position 
input angle is given by 
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(Eq. 7. 4. 13) 
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The response to an impulsive input is given by 



0 = 4 >. 



( 0 577 %) 



3 t 2 -0. 577 to t\ 
o 
e 



(Eq. 7. 4. 14) 



The response is seen to be well behaved at small gimbal 
angles; however, consider the location of the poles of the closed 
loop position control mode as the gimbal angle becomes larger. 
Figure 7. 4. 3 shows the migration of the poles for the value of 
K x chosen in equation 7. 4. 3. This figure indicates that at gimbal 
angles greater than 30 degrees the system becomes very oscilla- 
tory, and at approximately 84 degrees the system becomes unstable. 
The effect of the initial choice in open loop damping ratio, £ , 
is shown in Figure 7. 4. 4 which is plotted for the same K . A 
smaller would of course raise the curve of Figure 7. 4. 4 so 
that in the end a singular point at 90 degrees would result. 
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Figure 7.4.2 Steady state roll rate for constant 
sensor error signal. of 0.3 is an example only. 
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Figure 7.4.3 Location of Position Control Mode 
Closed Loop Poles for Large Gimbal Angles. 
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Figure 7.4 4 Plot of Boundary between Stable and 
Unstable Operation in Position Control Mode. 
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7. 5 Adaptive Mode 



It has been stated earlier that the four gyro system has 
the adaptive feature that permits continuous three axis control in 
the event that one of the four gyros is lost. This mode can be 
illustrated by assuming the loss of the number four gyro control- 
ler and examining the resulting spacecraft response. At the out- 
set, let it be assumed that the controller fails such that it loses 
its angular momentum at a steady rate over time, t^. Assuming 
no interaxial coupling this gives a moment disturbance about the 
roll axis as follows. 

—■ = M (Eq. 7. 5. 1) 

dt x ^ 

The moment disturbance is seen to be a step of H/t^ which lasts 
for t^ seconds, and this gives a total torque impulse of H. 

The spacecraft response to such a disturbance can be ob- 
tained from the equations of the previous section; however to 
illustrate the response for the adaptive mode the problem was 
set up on an analog computer and the results are shown in Figure 
7. 5, 1 and 7. 5. 2 for the case where a sensor is providing attitude 
errors. The zero input mode also operates to provide rate 
stabilization upon the loss of any controller. 
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Figure 7.5,1 Spacecraft Roll Rate and Roll Angle with Gimbal Rate 

and Angle for the Adaptive Mode, Data is from Analog Computer, 

I = 10 6 lb-ft-sec 2 , J = 10 lb-ft-sec 2 , t = 1000 sec. H = 10 4 
x 1 

lb-ft-sec, = 0.866. See Figure 7.5.2 for Phase Trajectory. 
Complete Run includes Initial Transient and Recovery Transient. 
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Figure 7.5.2 Phase Trajectory for Adaptive Mode 
I -10 6 lb-ft-sec 2 , J = 10 lb-ft-sec 2 , H = 10 4 lb-ft-sec. 
= 0. 866, See Figure 7.5.1. 
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7. 6 Summary 

This chapter has evaluated the response of the system 
to input disturbances and input commands for the zero input mode, 
the rate control mode, and the position control mode. The re- 
sults are shown in the illustration of the location of the poles of 
the various modes shown in Figure 7. 6. 1. This plot shows that 
an initial choice of the damping ratio of the zero input mode, £ , 
of approximately 0. 866 gives satisfactory performance for the 
rate control mode and the position control mode; therefore, no 
adaptive type change in the damping of the gimbal or tandem com- 
pensation is required to satisfy all three modes of operation. 

The effects of increase in gimbal angle is seen to make the zero 
input mode and the rate control mode less oscillatory, whereas 
the position control mode becomes more oscillatory. For a fixed 
value of gain, K there exists a gimbal angle for the position 
control mode which makes the system unstable. 

The system operating in the position control mode has 
inherent adaptive characteristics in that the failure of any one of 
the gyro controllers will automatically be compensated for by the 
repositioning of the opposite pair of controllers. 

The examples have been given for the roll equations. 
The pitch equations are almost identical except for the sign of 
the gain K^. The general solution for the yaw equations are 

identical except that the functions F. and F 0 have been defined 

IT o-y 

for a pair of controllers whereas is defined for a single con- 
troller. Also, the yaw controller uses all four controllers, and 
therefore faster response is available, if desired. 
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Figure 7,6,1 Illustration of the location of the poles of the Zero 
Input Mode (<^> =0.866), the Rate Control Mode (m=l), and 
the Position Control Mode, The poles are shown for gimbal angles 
of zero and 30 degrees to show effects of change in gimbal angle, 
The open loop poles of the position control mode are the closed 
loop poles of the zero input mode. 
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CHAPTER 8 



CONCLUSIONS AND RECOMMENDATIONS 

8. 1 Summary of the General Concepts o f the Thesis 

The thesis is concerned with the problem of attitude 
control of a spacecraft engaged in an extended mission. The 
foremost application of such a spacecraft is a manned explora- 
tion of the planet Mars. The factors considered in choosing a 
suitable attitude control system are 

Maximum reliability. 

Minimum ejection of mass. 

Minimum average power. 

Minimum system weight, and 

Minimum peak power. 

Attitude control systems which do not expend fuel mass 
but derive control moments by a time rate of change in angular 
momentum of a mass that remains within the spacecraft are 
called momentum exchange type attitude control systems and con- 
sume only power in their operation. In a space environment, 
power is more available and may actually be re- supplied from the 
sun, whereas expulsion control fuel is limited to that initially 
loaded aboard the spacecraft. Based on this, the main theme of 
the thesis is the study of momentum exchange type attitude con- 
trol systems. 

The thesis proceeds to accomplish five objectives set 
forth in section 1.2. In short, these objectives are 

derive the equations of motion, 

determine the torque disturbances. 



213 



select a specific control system, 

determine the response of the spacecraft, and 
present the conclusions. 

The attitude control system chosen to best satisfy the 
spacecraft requirements is a configuration of four controllers 
which operate to apply gyroscopic torques to the spacecraft. 

Each controller may be described as a gyroscopic type device 
having two degrees-of-freedom. The four controllers are ar- 
ranged in two pairs with each pair operating back-to-back. One 
pair of the controllers is actuated symmetrically to generate 
torques applied to roll the vehicle, whereas, the other pair, 
mounted normal to the first pair, provides torques for pitch con- 
trol. All four of the controllers provide torques for yaw control. 
This configuration of four controllers has a redundancy in its 
ability to effect a change in angular momentum along each of the 
three spacecraft control axes such that the system provides un- 
interrupted control upon the complete loss of angular momentum 
of any one controller. 

The four controller attitude control system is operated 
in three modes: 

Zero Input Mode 
Rate Control Mode 
Position Control Mode. 

The zero input mode enables the spacecraft to be rate 
stabilized in the absence of input commands. The rate control 
mode provides a means of changing the attitude of the vehicle at 
the maximum rate capability of the controllers. The position 
control mode is provided for alignment of the spacecraft with res- 
pect to a reference line of sight. Satisfactory stability charac- 
teristics were obtained for each mode. 
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8. 2 Summ ary of Chapter Con clusions 

Each chapter contains a brief summary, the highlights 
of which follow. 



Chapter 1: This chapter introduces the problem of attitude con- 

trol. In the practical problem the control system 
engineer must have knowledge of the physical charac- 
teristics of the spacecraft, the controllability re- 
quirements, and the stability requirements. 

Chapter 2: Chapter 2 contains the equations of motion of the 

spacecraft. The torques generated by the control 
system depend on the attitude of the controller 
rotor with respect to the spacecraft and the rates 
and accelerations of the controller rotor with res- 
pect to inertial space. The equations are presented 
in a manner which facilitates the evaluation of each 
of these variables. 



Chapter 3: One or more controllers are combined to form 

various attitude control systems. Control logic is 
required to give non-interacting control for the 
roll, pitch, and yaw vehicle attitude variables. 

A system of compensation is devised to minimize 
cross-coupling. 

Chapter 4: This chapter seeks to determine the torque dis- 

turbances acting on a spacecraft. It is determined 
that the external torque disturbances can be mini- 
mized by design and operation of the spacecraft. 
The torques which originate within the spacecraft 
from masses which do not leave the system have 
a zero mean value and their effects can be com- 
pensated for by using a momentum exchange type 
attitude control system. 



Chapter 5: 



Chapter 6: 



Chapter 7: 



This chapter compared the various systems and 
determined that the four gyro controller best ful- 
filled the control system requirements. The six 
gyro control system is the second choice. Inertia 
reaction wheel systems have very low power effi- 
ciency, and they lack the inherent stabilizing char- 
acteristics of the gyro controllers. Mass expulsion 
systems must be provided to some extent to de- 
saturate the momentum exchange system, but if 
used continuously, large amounts of fuel will be re- 
quired unless a limit cycle of the order of 15 min- 
utes is provided. 

The control loop was closed using the four gyro 
controller. A general solution of the equations re- 
quires machine computation because of the non- 
linearity and complexity of the equations. Three 
modes of operation are considered. Zero Input 
Mode, Rate Control Mode, and Position Control 
Mode. 

The response of the spacecraft is found to be satis- 
factory for all modes of operation using a fixed 
value of gimbal damping. The zero input mode and 
the rate control mode become less oscillatory at 
larger gimbal angles, whereas the position control 
mode becomes more oscillatory. The error signal 
from the sensor must be limited to avoid a possible 
unstable condition which drives the gimbal to their 
full 90 degree position. 
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8. 3 Summary of Characteristics of Gyro Type Controllers 

The following conclusions are presented concerning 
gyro type controllers used to provide attitude control of space- 
craft. 

a. The use of gyro controllers operating back-to-back 
eliminates large cross control moments typical of single gyro 
controllers, and this becomes more important as the gimbal 
angle is increased. 

b. Gyro controllers should be operated without an ap- 
preciable amount of saturation in order to minimize gyroscopic 
cross coupling moments. This statement is true about any mo- 
mentum exchange system, and suggests that zero angular mo- 
mentum type systems are preferred to those which do not have 
zero angular momentum in their initial configuration. 

c. Except for the unstable, full 90 degree position of the 
gyro controllers where no control moment is required, gyro 
controllers require continuous control moments to provide an 
attitude rate to the spacecraft. 

d. During the angular acceleration of the spacecraft to 
an attitude rate the torque multiplication for the gyro torquer is 
high. The ideal torque gain for a pair of controllers operating 
open loop is equal to 

1 

*o 

e. Gyro controllers provide inherent rate stabilization. 

f. Since the gyro gimbal is free to change its attitude 
with respect to the spacecraft, generally a control system using 
several gyro controllers will have inherent adaptive characteris- 
tics, in that, the failure of a single gyro unit will be compensated 
for by the repositioning of one or more of the other gyros. 
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g. In the absence of external torques the spacecraft as- 
sumes a particular attitude rate corresponding to a controlled 
gimbal position as given by the following equation for roll rate. 



p sin 7 .. 
*max 1 



n. Gyro controllers are more efficient than wheel con- 
trollers since they can effect a change in angular momentum of 
the spacecraft without a change in their kinetic energy, whereas 
wheel controllers require a kinetic energy change to effect an 
angular momentum change which results in a wheel efficiency 
ratio of the order of the ratio of the wheel moment of inertia to 
the spacecraft moment of inertia. 

i. The provision of a position control loop for the gimbal 
angle, as is done in the rate control mode, causes a deterioration 
in the ability of the gyro to provide inherent rate stabilization. 

j. The ratio of moment of inertia to damping coefficient 
for the gyro controller should be greater than that typically used 
for integrating gyros of inertial navigation systems, 

k. The power requirement for a gyro controller is ap- 
proximately constant as compared to an inertia wheel which has 
severe peak power requirements. 

l. Gyro controllers can provide fast, accurate, and well 
damped control for manned spacecraft. 
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8. 4 Recommendations for Further Study 

As an extension to the thesis it is considered desirable 
to continue the study in the following areas. 

a. Optimization of the design of gyro controllers. 

b Study of controllers which employ fluids to obtain a 
change in angular momentum, 

c. By the use of machine computation, determine the 
effects of the cross coupling moments acting on the spacecraft. 

d Determine the feasibility of using existing angular 
moment such as power turbines in a spacecraft to achieve atti- 
tude control. 

e. Investigate devices which can store energy as well as 
deliver energy for use with inertia reaction wheel control sys- 
tems . 
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APPENDIX A 



SYMBOLS AND MATRIX NOTATION 



A. 1 Matrix Notati on 

The treatment of an analytical problem that uses several 
sets of orthogonal Cartesian coordinates is simplified if matrix 
methods are employed so that the coordinate frames are unam- 
bigously inferred by the matrix expressions. In reference 1 there 
is contained a particularly good explanation of matrix notation as 
applied to control system problems. Thus for detailed explanation 
the reader is referred to that reference and only a brief discussion 
of the notation is presented herein. 

A.2 Coordinate Transformation Matrices 

The symbol Q is defined as an orthogonal transformation 
between two Cartesian coordinate frames, and it is a 3 X 3 square 
matrix. Each Q shall contain a double subscript indicating the 
pertinent coordinate reference frames in a "to-from" sequence 
from left to right. For example, Qj^ is the coordinate transfor- 
mation which is post multiplied by a vector in the A frame to 
transform the vector to the I frame. Q^j then is the coordinate 
transformation that transforms a vector from the I frame to the 
A frame, and since all Q transformations are orthogonal trans- 
formations the operation of inverting a Q matrix is simply that of 
taking the transpose of the matrix. 




(Eq. A. 2. 1) 



Choice of the "to-from" sequence in subscripts facilitates 
multiple transformations so that the following is true. 

Q IR = Q IV Q VA q ar (Eq ‘ A ' 2 ‘ 2) 
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It is sometimes necessary to represent a coordinate frame 
by more than one letter, as for example the gimbal frame is re- 
ferred to as GIM. When there is chance of ambiguity in a subscript 
a comma is used to separate the two frames. Thus the coordinate 
transformation Qjjyj represents the coordinate transformation 
to the A frame from the GIM frame. The coordinate frames used 
in this thesis are contained in Appendix C. 

A. 3 Relative Velocities 

A velocity is a vector quantity and is expressed in matrix 
form by a column vector. Thus Wj R means the 3X1 matrix rep- 
resenting the angular velocity of frame E relative to frame I, ex- 
pressed in frame E. 



In matrix equations it is necessary to have a means of 
expressing the operation of a vector cross product. Thus consider 
a vector 



W^ R = li+mj + nk 



(Eq. A. 3.1) 



and a vector 



H r = Ai + Bj + Ck 



(Eq. A. 32) 



The cross product W TD X H 

IK 

(Refer to ref 29 
p 190) 




1 

A 



j 

m 

B 



k 



n 

C 



(Eq. A. 33) 



W IR X H r = (m C- nB) i > (nA-lC)j + (lB-mA)k 

(Eq. A. 34) 



In matrix form the result of the cross product operation can be 
represented as a column matrix such that. 



w 

IR 



h rJ 



R 



m C - n B 
nA - 1C 
IB - m A 
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(Eq. A. 3.5) 



Then upon separation of the variables associated with 



Wjp, and 



d 



R 



respectively gives: 




(Eq. A. 3.6) 



Therefore we define from any column matrix such as 

1 



W IR = 



the 3X3 antisymmetric matrix 



m 



n 



(Eq. A. 3.7) 



W IR ^ 5 



0 

n 

•m 



-n m 

0 -1 

1 0 



(Eq. A. 3.8) 



Reference 1 derives this definition in a more rigorous 
manner from the fundamental operation of differentiating a matrix 
equation. Suppose 

W IR = SRA W IA + W AR (E <>- A ' 39> 

Differentiating with respect to time 

Wir = «RA ™IA + <5rA W !A + War ®q. A. 3.10) 



Since Qp>^ ^AR = ^ a ma t r i x Qj^a can ^ ac ^ ore< ^ from 
equation A. 310. 



W IR = QRA 



W IA + Q AR Q RA W IA 



+ w ar (Eq. A. 811) 
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and by differentiating each of the elements of Q . and premulti- 
plying the result by we find that 



of the form 



Qj>a * s a 3 x 3 antisymmetric matrix 



= as defined in equation A. 3.8. 
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A. 4 Numbering of Controller Elements 

When control systems are considered that contain many 
elements some system must be devised to keep track of the 
various elements. For the lack of a better arrangement the con- 
trollers have simply been numbered as the problems have been 
solved. Therefore, the following is a list of the various numbers 
assigned to the controllers (Gyro). 

1 Gyro with Spin Reference Axis along y axis 

2 Gyro with Spin Reference Axis along -y axis 

3 Gyro with Spin Reference Axis along x axis 

4 Gyro with Spin Reference Axis along -x axis 

5 Gyro with Spin Reference Axis along z axis 

6 Gyro with Spin Reference Axis along -z axis 

7 Gyro with Spin Reference Axis in x-y plane 

and rotated 120° from Gyro 1 about z axis. 

8 Gyro with Spin Reference Axis in x-y plane 
and rotated -120° from Gyro 1 about z axis. 

9 Gyro with Spin Reference Axis in y-z plane and 
rotated -o’ degrees from -y axis about x axis 

10 Gyro with Spin Reference Axis in y-z plane and 
rotated I -a degrees from -y axis about x axis 

A particular control system may require single -degree - 

of-freedom controllers, two-degree-of-freedom controllers, or 

possibly three-degree-of-freedom controllers in which the angular 

speed of the rotor may be varied in addition to the two degrees- 

(33) 

of-freedom in precession of the spin axis of the controller 
Therefore the controllers used for a particular system require 
further specification to fully define their configuration. Accord- 
ingly in Appendix G which presents specific control systems the- 
angles required to align the case of the controller shown in 
Figure B. 9 of Appendix B are given for each system. See also 
Appendix C. 3 for coordinate transformations for particular con- 
trollers. 
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Figure A. 43 illustrating general positions of 
Controllers 9 and 10. 



A. 5 List of Symbols 



a Angle used to Define Gimbal Position in Case. See Figure 

B. 10. 

y Angle used to Define Case Position Relative to Vehicle. 

See Figure B. 9. Also used as a General Control System 
Variable. 

j3 Angle used to Define Rotor Position Relative to Gimbal. 

See Figure B. 11. 

■q , n , & n Represent Primary Control Variables. Also used 
x y z 

as angle without subscript in Figure B. 3. 

fp Represents Vernal Equinox. 

&i Angular Rotations of Heliocentric Orbital Plane 

Reference Frame relative to Heliocentric Inertial 
Reference Frame. 

<f>, 6 and Angular Rotations of vehicle-centered Vehicle 

Reference frame with respect to an Inertial 
Reference Frame. 

^K'^K’k^KP Refers to Angular Rotations of Geocentric Orbital 
Position Reference Frame relative to the Geo- 
centric Orbital Plane Reference Frame. See 
Figure B. 13. 

LAT Latitude. See Figure B. 14. 

LON Longitude. See Figure B. 14 

p and A Angular Rotations used in Defining the Geocentric 

s s 

Solar Reference Frame. See Figure B 15. 

A Represents an Angle when provided with a subscript. When 

used as a subscript the symbol refers to the Vehicle 
Principal Axis Reference Frame. See Figure B. 8. 

E Energy 
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H 

M 

P 

Q 

R 



p> q. 



w 



LA. 



I 

J 

n 

• 




S 

c 



Angular Momentum 

Moment 

Power 

Orthogonal Coordinate Transformation Matrix. (See Ap- 
pendix C) 

Resistance, ohms 

and r Vehicle Attitude Rates in Roll, Pitch, and Yaw 

Angular Velocity of Coordinate Frame A relative to 
Coordinate Frame I expressed as a Column Matrix 

Moment of Inertia of Spacecraft 

Moment of Inertia of Components of Control System 

Refers to number of columns or rows of a matrix 

A dot over a symbol indicates a differentiation with re- 
spect to time 

Represents the forming of a Matrix from a Column Vec- 
tor. See Section A. 3. 

Shorthand notation for the trigonometric sine. 

Shorthand notation for the trigonometric cosine. 



Subscripts 
T Total 

CS Control System 

ext External 

VF Refers to Masses in the Vehicle that are Rigidly Fixed 

to the Spacecraft 

VM Refers to Masses in the Vehicle that Move Relative to the 
Vehicle Exclusive of the Rotating Member of the Control 
System Controller. 
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N 



c 

g 

r 



Represents a summation symbol for N Controllers 
Armature 

Summation Index. Also used in defining the Heliocentric 
Orbital Plane Reference Frame. See Figure B. 2. 

Refers to the Case of the Controller 

Refers to the Gimbal of the Controller 

Refers to the Rotor of the Controller 



x, y, and z Generally refers to the x, y, and z direction of an 
Orthogonal Coordinate Reference Frame. 

w Wheel 

Z <[>, YN, and XV Refers to Angular Rotations of Vehicle-Cen- 
tered Vehicle Reference Frame relative to 
Vehicle-Centered Inertial Reference Frame. 
See Figure B. 7. 

ZV, 1; YY;andXA Refers to Angular Rotations of Vehicle-Cen- 
tered Principal Axis Frame relative to 
Vehicle-Centered Vehicle Reference Frame. 
See Figure B. 8. 

YU and X, GU Refers to Angular Rotations of Vehicle-Cen- 

tered Gyro Case Axis Reference Frame rela- 
tive to Vehicle-Centered Vehicle Reference 
Frame. See Figure B. 9. 

Z1 Refers to angular rotation of Earth. See Figure B. 12. 



See section B. 1 for Symbols used for Coordinate Refer- 
ence Frames. 
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A. 5 Glossary of Terms 



Active Control 

Active control of a spacecraft is defined as the operation 
of controlling the spacecraft with a torque producing con- 
trol system operating with error sensors in a closed loop. 
The control system always 'consumes energy. 

Compensation 

A signal provided to a control loop which is proportional 
to a signal of another loop, and is used to minimize the 
effects of an unwanted coupling between these two loops. 
Three types of compensation are used in this report: 
gyroscopic coupling, cross control coupling, and space- 
craft inertia cross coupling. 

Control Logic Matrix 

The control logic matrix is defined as an n x 3 matrix 
which pre-multiplies the primary control variables to 
define individual signals to the n-degree-of- freedom con- 
troller. 

Control System Coupling Matrix 

The control system coupling matrix is a 3 x 3 matrix 
which operates on the vehicle rate variables resulting 
from the combination of one or more controllers. If the 
control system coupling matrix contains no diagonal terms 
then it is also the gyroscopic coupling matrix. 

Control System Input Matrix 

A 3 x n matrix which results from any arrangement of 
terms of one or more controllers in which the control 
system input matrix operates on the control systepn input 
variables. 
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Control System Input Variables 



Any controller may have one or more degrees of freedom 
which may be considered the input variables of the con- 
troller, i. e. with wheels, the input variables are angular 
acceleration of the wheels. With pure gyro systems the 
input variables are precession rates. A combination of 
two or more controllers give many degrees of freedom 
which represent the control system input variables. 

Controller 

A momentum exchange device which is capable of applying 
control moments to the spacecraft. A controller in this 
report is considered to be a rotating rigid body. One or 
more controllers are combined to form a complete con- 
trol system. 

Gyroscopic Coupling Matrix 

Defined as the control system coupling matrix with all 
diagonal elements replaced by zeros. 

Guidance 

The guidance of a vehicle is defined as the operation of 
controlling the thrust vectors acting on a spacecraft, such 
that, a desired trajectory is followed. Guidance of a 
spacecraft is required during thrusting and possibly during 
re-entry if the aerodynamic lift of the vehicle can be con- 
trolled. 

Inertial Guidance Measurement Unit 

The inertial guidance measurement unit is defined as an 
assemblage of instrumentation to determine the specific 
force vectors acting on the spacecraft suitable for per- 
forming the guidance and navigation function. 
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Momentum Exchange 



A momentum exchange control system is defined as a de- 
vice which is capable of providing control torques to a 
spacecraft by a time rate of change in angular momentum 
within the device. Consequently, the total mass and angu- 
lar momentum of the combined vehicle and control system 
remains constant. 

Momentum Transfe r 

A momentum transfer control system is defined as a de- 
vice which is capable of providing control torques to. a 
spacecraft by ejecting mass from the spacecraft and 
creating a force-impulse normal to a lever arm directed 
to the center of gravity of the spacecraft. 

Navigatio n 

Navigation of a spacecraft is defined as the operation of 
determination of the position and velocity (or position and 
velocity deviations), and computation of the guidance com- 
mands necessary to arrive at the desired destination. 

Re st Po int 

A spacecraft rest point is. the attitude which results in 
zero applied torque, and at which attitude the spacecraft 
exhibits stable static Stability characteristics. 

Passive Contro l 

Passive control of a spacecraft is defined as the operation 
of controlling the spac'ecraft purely, by means of an exist- 
ing stabilizing torque which acts on the spacecraft. At- 
titude control is achieved by dissipating energy and the 
system provides inherent sensing. 

Pointing Accuracy 

The pointing accuracy is defined as the maximum deviation 
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of the output of a closed loop with reference to a desired 
output when no disturbances are admitted to the loop. 
Pointing accuracy is a measure of the stability of a 
spacecraft controlled by an active control system. 

Primary Control Matrix 

The primary control matrix is a result of pre- multiplying 
the control logic matrix by the control system input ma- 
trix. The primary control matrix is then a 3 X 3 matrix 
and is diagonal or nearly diagonal for the range of control 
system input variables chosen to define the control logic 
matrix. 

Primary Control Variables 

Three variables which give non-interacting (or nearly 
non-interacting) control in roll, pitch, and yaw. 

Rest Position 

For a spacecraft in an environment in any time invariant 
situation, there exists a spacial orientation at which the 
external moments acting on the spacecraft exactly cancel 
and no moment acts on the vehicle. 

Semi-Passive Control 

Semi-passive control of a spacecraft is defined as passive 
control systems which increase their damping properties 
by introducing gyroscopic effects, or those systems which 
possess no static stability but have strong damping charac- 
teristics such as spin stabilized vehicles. 

Saturate (and desaturate ) 

Any momentum exchange type control system will have 
an initial angular momentum disposition and a maximum 
angular momentum storage capability for a particular 
control axis. Saturation is defined as the percent change 
in the angular momentum of a particular axis when the 
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vehicle is non-rotating with respect to the reference frame. 
Desaturation is defined as providing an external moment 
to the spacecraft which tends to return the momentum 
exchange control system to its initial disposition. 

Vehicle Attitude Rate Variables 

The vehicle attitude rate variables are the roll, pitch, 
and yaw rates of the vehicle with respect to inertial 
space. 



234 



APPENDIX B 



COORDINATE REFERENCE FRAMES 



B. 1 


Summary of Coordinate Reference Frames Defined 


Figure 


Symbol 


Title 


B. 1 


II 


Heliocentric Inertial Reference Frame 


B. 2 


H 


Heliocentric Orbital Plane Reference Frame 


B. 3 


B 


Heliocentric Orbital Position Reference 






Frames 


B. 4 


<£, 0,0 


Vehicle-Centered Solar Orbital Reference 






Frame 


B. 5 


I 


Vehicle-Centered Inertial Reference Frame 


B. 6 


VR 


Vehicle-Centered Velocity Reference Frame 


B. 7 


V 


Vehicle-Centered Vehicle Reference Frame 


B. 8 


A 


Vehicle-Centered Principal Axis Frame 


B. 9 


GU 


Vehicle-Centered Gyro Case Axis Reference 






Frame 


B. 10 


GIM 


Vehicle- Centered Gyro Gimbal Axis Reference 






Frame 


B. 11 


R 


Vehicle-Centered Gyro Rotor Axis Reference 






Frame 


B. 12 


III 


Geocentric Inertial Non-rotating Reference 
Frame, and 




E 


Geocentric Earth Reference Frame 


B. 13 


K 


Geocentric Orbital Plane Reference Frame, 




P 


and Geocentric Orbital Position Reference 






Frame 
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B. 14 


G 


Geocentric Longitude- Latitude Grid Reference 
Frame 


B. 15 


S 


Geocentric Solar Reference Frame 


B. 16 


o 


Vehicle-Centered Planet Orbital Reference 
Frame 



An attempt has been made to keep the defined reference 
frames identical with those given by Ogletree in reference 21. 
For reference frames centered at a planet of the solar system 
other than earth, it is considered that Figures B. 12 and B. 16 
can be applied with a suitable subscript denoting the planet con- 
cerned. 
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HELIOCENTRIC INERTIAL REFERENCE FRAME II 







r 

VERNAL EQUINOX 



FIGURE 8. 1 An illustration of an iner tia! fixed frame assuming the sun as a fixed point in space 



HELIOCENTRIC ORBITAL PLANE REFERENCE 
FRAME H 

AXIS X H PASSES THROUGH PERIHELION OF SPACECRAFT ORBIT 




Q 



H,n : 



(CwCfl-SwSiiClKCwSil + SwC&ClMSwSI) 
(-SwCfl - CwSiiCiK-SwSii + CwCflC I HSICoi) 
SiSii -SlCil Cl 



FIGURE B. 2 An Illustration of on Orbital Plane Reference Frame defined as Three Rotations from the 
Heliocentric Inertial Reference Frame. 
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FIGURE B. 3 An illustration of a Heliocentric Orbitol Position Reference Frome defined by a Single 
Rotation from the Heliocentric Orbital Plane Reference Frame. 
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VEHICLE-CENTERED SOLAR ORBITAL REFERENCE 
FRAME (THREE CASES) 






-IOO 




0 1 0 




0+1 0 


°<*> 8 s 


0+1 0 
0 0-1 


°e B 5 


1 0 0 
0 OH 


it 
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O 


To 
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FIGURE B.4 An Illustration of a Vehicle-Centered Orbital Reference Frame asreloted to the Heliocentric 
Orbital Position Reference Frame. 
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VEHICLE-CENTERED VEHICLE REFERENCE FRAME 3 Z 
SHOWN RELATIVE TO VEHICLE-CENTERED INERTIAL 
REFERENCE FRAME I 




1. Frame I is non-rotating with respect to fixed stars, i.e. with respect to 

FRAME 3T. 

2. Directions of Xj ,Yj ,Zj are defined to coincide with X 0t Y 0 ,Z 0 at time 

ZERO FOR ANY PARTICULAR SET OF INITIAL CONDITIONS. 

3. Order of rotations to place frame i in coincidence with frame zl. 



A. 


Rotate about Axis Z x 


through Angle if/. 


B. 


Rotate a bout Axis Y v 


through Angle 9. 


C. 


Rotate about Axis X v 


through Angle <f>. 



FIGURE 0-5 An illustration of relation between the Vehicle-Centered Inertial Reference Frame and 
the Vehicle -Centered Vehicle Reference Frame 
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VEHICLE -CENTERED VELOCITY-REFERENCE 
FRAME VR 




WvR ,(£ = 




QvR,<£ 



+CA Z0 +SA Z0 0 
-SA Z0 4- CA Z0 0 
0 0 I 



FIGURE B.6 An illustration of a Vehicle -Centered Velocity- Reference Frame. Angle A Z q is 
rotation about Z to place Framel in coincidence with Frame VR. 
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VEHICLE-CENTEREO VEHICLE REFERENCE FRAME 3ZL 



DIRECTION OF 




Order of rotations to place Frame < f > in coincidence with Frame 3ZL . 

1. Rotate about Axis through Angle A Z( ^ 

2. Rotate about Axis Y n through Angle A Y n* 

3. Rotate about Axis X v through Angle A xv . 

See Appendix C for Q v< £. 



FIGURE B . 7 An illustration of a Vehicle - Centered Vehicle Reference Frame as defined by three 
rotations from the Vehicle -Centered Orbital Reference Frame. 
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FIGURE B.8 An Illustration of a Vehicle-Centered Principal Axis Frame as defined by three rotations 
from the Vehicle* Centered Vehicle Reference Frame. 
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VEHICLE-CENTERED GYRO CASE AXIS REFERENCE 
FRAME GU 




Order of rotations to place Framed in coincidence with Frame GU. 

1. Rotate about Axis Z v through Angle 

2. Rotate about Axis Yy through Angle A yu . 

3. Rotate about xis X 6lj through Angle A x gu * 

See Appendix C for Q gu 



W V,GU 


1 

o o 
1 


J Cgu = 


Jcx 

0 


0 

J CY 


O O 
1 




r 




0 


0 


J cz 



FIGURE B.9 An illustration of a Vehicle-Centered Gyro Case Axis Reference Frame os defined by three 
rotations from the Vehicle-Centered Vehicle Reference Frame. 
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VEHICLE-CENTERED GYRO GIMBAL AXIS REFERENCE 
FRAME GIM 




To place Frame GU in coincidence with Frame GIM. 
I. Rotate about Axis X gu through Angle a . 





a 




J,X o o' 


W GU,GIM = 


0 


J 9 . 


GIM = 


0 J gY 0 




0 




0 0 JgZ 




1 


0 0 




Qqim,gu = 


0 


4* cos a + sin a 






0 


-sin a +cosa 





FIGURE B.IO An Illustration of a Vehicle -Centered Gyro Glmbal Axis Reference Frame as defined by a 
single rotation from the Vehicle-Centered Gyro Cose Axis Reference Frame. 
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VEHICLE-CENTERED GYRO ROTOR AXIS REFERENCE 
FRAME R 




To place Frame GIM in coincidence with Frame R 
I. Rotate about Axis Y 6im through Angle /3. 





0 




Jrx 


0 0 


W GIM,R = 


9 

P 


JrR = 


0 


J r v 0 




0 




0 


0 Jrz 



Qr.gim 



cos/3 O -sin/3 

0 I 0 
+sin/3 0 cos/3 



FIGURE B.ll An illustration of a Vehicle- Centered Gyro Rotor Axis Reference Frame as defined by a 
single rotation from the Vehicle* Centered Gyro Gimbal Axis Reference Frame. 
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GEOCENTRIC INERTIAL NON-ROTATING REFERENCE 
FRAME IE AND GEOCENTRIC EARTH REFERENCE 



FRAME E 



| z m> z E 






CA zr 


SAzi 


0 




0 




0 


°E,zr 5 


-SA ZI 

o 


CA zr 

n 


0 

| 


w hl,e 1 


0 

2 TT RAD 


z 


0 

"7 'tftone Ipk ”5 RAO 




\J 


V 


1 




86,164 SEC 




f.cyciio xiu 5 £Q 



FIGURE 6*12 An Illustration of the relations between a Geocentric Inertial Non-Rotating Reference 
Frame ond a Geocentric Earth Reference Frame which rotates at Eorth Rote. 
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GEOCENTRIC ORBITAL PLANE REFERENCE FRAME K 
AND GEOCENTRIC ORBITAL POSITION REFERENCE 




X K : -(Wxk + W ie XT-To ) + X*J 0 
T = INITIAL VALUE OF TIME 

X K J 0 = INITIAL LONGITUDE OF ASCENDING NODE OF ORBIT 



FIGURE B.I3 An illustration of the relations between a Geocentric Orbital Plane anda Geocentric 
Orbital Position Reference Frame 
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GEOCENTRIC LONGITUDE-LATITUDE GRID 
REFERENCE FRAME G 




Y IS IN DIRECTION OF 
W 1E * R P 



ClatClON 


ClatSlon 


Slat 


" S LON 


C L0N 


0 


"SlatClon 


“ S LAT C LON 


Clat 



FIGURE B.I4 An Illustration of a Geocentric Longitude-Latitude Grid Reference Frame. 
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GEOCENTRIC SOLAR REFERENCE FRAME S 




fj-s ^ s Jo + W £S (T-T 0 ) 

/is | 0 - INITIAL VALUE OF /i s 



FIGURE B.I5 An Illustration of a Geocentric Solar Reference Frame. 



251 



VEHICLE -CENTERED PLANET ORBITAL REFERENCE 
FRAME 0 




w k 0 Hw kp 

$o3y p 

Yq£ 3 5 p Q 0 p ■ 

A A 

z 0 HXp 



o i o 
o o-i 
-10 0 



FIGURE B.I6 An illustration of a Vehicle-Centered Orbital Reference Frame. 
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APPENDIX C 



COORDINATE TRANSFORMATIONS 



C. 1 Coordinate Transformations Between Reference Frames 
Relating to Interplanetary Space Analysis 

(Refer to Appendix B for definitions of coordinate frames) 



Q 



H, II 



Q 



Q 



0 B 



(CwCfi - SuSfiCi) (CuSfi + SuCnCi) ScoSi 

(-SuCfi-Cusnci) (-Susn + Cucnci) siCw 



BH = 



S iSQ 



Cn 


Sn 


0 


-S n 


C n 
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0 


0 


1 


-1 


0 


0 


0 


1 


0 


0 


0 


-1 






— ’ 



-sicn 



Ci 



(Eq. C. 1. 1) 



(Eq. C. 1. 2) 



See Figure B. 4 for 

J 0B 

(Eq. C. 1. 3) 



«0B and 



Q 



VR, 0 



CA zo 


SA zo 
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“ SA zo 


CA zo 
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(Eq. C. 1. 4) 
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(Eq. C. 1.7) 
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TRANSFORMATIONS FOR SPECIAL CASE 
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C. 2 Coordinate Transformations Betwe en Reference Frames 
Relating to Planetary Space Analysis 
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NOTE: 



C. 3 Coordinate Transformations Relating to the Definition of 
Particular Controllers 
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SINGLE -DEGREE- OF- FREEDOM CONTROLLERS 



Q 



GUI, V 



1 

0 

0 



0 

1 

0 



0 

0 

1 




Q GU2, V 



-1 0 

0 -1 
0 0 




(Eq. C. 3. 8) 



263 



1 



Ygu 



Q GU3, V " 



0 0 

1 0 0 

0 1 o 




Vv 

£* (Eq. C. 3. 9) 



Q, 



GU 4, V 



■1 

0 



0 

1 



Q, 



GU 5, V 



0 

1 



0 

0 



Q, 



GU 6, V 



0 

0 

- 1 



+1 

0 

0 





264 



0 



0 



Q GU8, V 



Q GU9, V 



Q GU10, V 



+ S 7 + C 7 
o 'o 

-C 7 +S 7 

o 'o 

-1 0 

0 -Ca 

I 0 -So? 

L 




(Eq. C. 3.15) 



-1 0 

0 -Ca 

o 

0 +Scv 

o 




265 



APPENDIX D 



D. 1 Relative Velocities 
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r (Eq. D. 1. 2) 

r^ = radius of perihelion of spacecraft orbit 
= velocity of spacecraft at perihelion 
r = radius of spacecraft from sun 

(Eq. D. 1. 3) 
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D. 2 



Approximate Relative Velocities 



Valid when the frames nearly coincide or when angular 
velocity is purely about a single axis. The subscript A may be 
substituted for V for the case where Qy ^ = I. 
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for the particular problem where 
|3 » a and j3 » 7 which assumes the 
rotor spin as large. 



D. 3 



Relative Velocities for Controller 



For the special case where the controller is in position 
1 , i. e. A^. Qy = 0 , Ayy = 0 , and 7 is a control variable , 
the following relative velocities apply. 
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D. 4 



Time Rate of Change of Angular Velocities 

* • • • • • 

For the case where p » tt, y, p, q and r. Controller 1 
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APPENDIX E 



EXACT MOMENT EQUATIONS FOR 
CASE, GIMBAL AND ROTOR TERMS 

In accordance with the definition of the GU, the GIM, 
and the R coordinate reference frames of Appendix B, the case 
containing a gimbaled rotor is chosen so that the spin reference 
axis of the rotor is aligned along thb positive y axis of the V frame. 
This gyro position is also called gyro 1. It is assumed that the 
V frame (vehicle- centered vehicle reference frame) and the 
A frame (vehicle- centered principle axis frame) coincide. 

Equation E. 3. 1 is applicable to problems concerning 
inertia reaction wheels, but it is not particularly suited for 
wheels. Accordingly, section F. 2 of Appendix F is a derivation 
of an equation specifically for inertia reaction wheels. 
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APPENDIX E 



EXACT EQUATIONS FOR CASE, GIMBAL, AND ROTOR TERMS 



E. 1 Case Terms 

From Equation 2.3.6 
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Substituting these matrices and performing the indicated 
operations gives for the contribution of the case terms the 
following moments. 
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(Eq. 



E. 2 



Gimbal Terms 



By an identical procedure as the previous section, the 
gimbal terms are found to be the following. 
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G imbal Terms 
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-^(pS7 " q C7) (pC7 + qS7 + a) 
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rpc<i>ce - es <t> 



The preceding is adequate when two of the three attitude 
angles are small, otherwise we must substitute 



p = i - i se 
q = 0 c 0 - is <j)Ce 
r = i/jC<j)C 0 - 0S<t> 
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E. 3 



Rotor Terms 



By an identical procedure as in section E. 1 the exact 
rotor terms for a rotor of a controller in position 1 are found 
to be the following. 

Define 

J i dimetrical moment of inertia of rotor 
rd 

J = polar moment of inertia of rotor 
rp ^ 

The following equations are adequate when two of the 
three attitude angles are small. Otherwise we must substitute 

p = (j) - ipS6 

q = 9 C<t> - is<t>C6 

r = 0C0C0 - 6S<t> 
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Rotor Terms 
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APPENDIX F 



APPROXIMATE EQUATIONS FOR CONTROLLERS 



F. 1 Arbitrary Controller 

Using Equation 2. 3. 8 and assuming the angular momentum 

(33) 

of the rotor is predominately along the spin axis of the rotor, 
one can derive the approximate equation 



Ksj Q A, GIM W A, GIM*^ H rJ 



GIM 



+ W IA^°A, GIM H rJ 



GIM 

(Eq. F, 1. 1) 



where H 
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GIM 



Jp0 



(Eq. F. 1. 2) 



Now W 



IA 
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q 

r 



from Appendix D 



(Eq. F. 1. 3) 



Q A, GIM °A, GU °GU, GIM 



(Eq. F. 1. 4) 



and can be found from the 
coordinate transformations 
of Appendix C. 1 



W A, GIM °GIM, GU W A, GU + W GU, GIM 

(Eq. F. 1. 5) 

and can likewise be found 
using matrices of Appendix C 
and Appendix D. 
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Note that equation F. 1. 1 is a sum of two parts. The first 
of these parts will form the primary control matrix and the 
second part will form the gyroscopic coupling matrix of the final 
equation in a particular system using one or more of the control- 
lers. 

After performing the indicated operations, the equation for 
a single controller with an arbitrary case position is given by the 
following equation. 
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It is not considered desirable to put the preceding equation 
in matrix form as there would be no simplification. However, 
upon use of the preceding equation in a particular controller 
position, the time rates of the case angles will usually vanish and 
considerable simplification then results. At this point it is de- 
sirable to put the remaining terms in matrix form to assist their 
summation with terms from other controllers with the same 
control system input variables. 
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F. 2 



Equations for Inertia Wheels 



Unfortunately equation F. 1. 6 cannot be used for the case 
of inertia reaction wheels because it has been assumed in that 
section that the angular velocity of the controller is a large value 
and is constant. Equation E. 3. 1 can be used for the X wheel and 
the Y wheel or any other wheel that lies in the x-y principle axis 
frame or for a pure z -aligned wheel. For the perfectly arbitrary 
inertia reaction wheel it may be desirable to list here a general 
expression for the moments. If we consider that the wheels are 
rigidly mounted in the spacecraft then equation 2.3.8 can be used 
and is written as follows. 
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Or, the same equation written with respect to the case frame is 
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Since the case holding the wheel is rigidly attached to the space- • 
craft W AGU ^= 0 
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(Eq. F. 2. 2) 



293 



The most difficult part of this derivation is that of finding 
H„ i . A considerable simplification of the algebra results 
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if we agree in this section to let 
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(Eq. F. 2. 3) 



Where the elements of equation F. 2. 3 are given by equation C. 1. 7 
for the case where the vehicle reference frame is identical with 
the principle axis frame. Then it is determined that the angular 
momentum with respect to the case frame is as follows. 
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Following through the indicated operations by substituting 
the above matrices into equation F. 2. 2 gives the following equa- 
tion for the moments contributed by the inertia reaction wheels 
of an attitude control system. 
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(Eq. F. 2. 7) 



Equation F. 2. 7 has been written in two parts because it is 
convenient to lump the second part with the vehicle since that part 
is not dependent on the control variable /3. The second part relates 
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to the moments contributed to the vehicle by the wheels acting as 
inert masses constrained to move with the vehicle. Therefore the 
equation is written as the following. 
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from equation F. 2. 3 and C. 1. 7 
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For a Y-wheel 
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F. 3 Equations for a Particular Three Degree-of- Freedom 
Controller 



Only one controller is considered wherein three degrees- 

of- freedom of the rotor are allowed. This controller is aligned 

with its spin reference axis along the positive x-axis of the 

vehicle-centered principal axis frame. Thus in equation E. 3. 1 

if we let 7 - 7-90° and J = 0, the following equation is obtained 
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(Eq. F. 3. 1) 
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F. 4 



Equations for Particular Two- Degree- of- Freedom 
Controllers 



In the synthesis of the moment equations for those control 
systems that are composed of several controllers of the'two- 
degree-of- freedom type, it is convenient to have the simplified 
moment equations of each controller. Thus the controllers de- 
fined by the following angles are presented. See Figure B. 9 and 
paragraphs A-4 and C-3. 



TABLE F. 4 

POSITION ANGLES FOR TDF CONTROLLERS 



TDF CONTROLLER 


CASE ANGLES 


A X, GU 


a yu 
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0 +y 
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0 
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± 180° + y 
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"o 
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- 90° + y 
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+ 90° + y 
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O' 


+ 120° + y 


8 


0 


0 


- 120° + y 



If the above reference angles are substituted into equation 
F. 1. 6 the following equations are obtained. Note that y is a case 
angle retained as a control input variable. 
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(Eq. F. 4. 3) 
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(Eq. F. 4 . 5 ) 
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F. 5 Equations for Particular Single- Degree- of- Freedom 

Controllers 

In the synthesis of the moment equations for those 
control systems that are composed of several controllers of the 
single-degree-of- freedom type, it is convenient to have avail- 
able the simplified equations of each controller. The following 
table of angles define the orientation of the case relative to the 
vehicle centered principal axis frame. See Figure B. 9. 

TABLE F. 5 



POSITION ANGLES FOR SDF CONTROLLERS 



SDF CONTROLLER 


• CASE ANGLES 
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SINGLE- DEGREE- OF- FREEDOM CONTROLLERS 
(Subscripts of Variables Omitted)* 



Controller 1 
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(Eq. F. 5. 1) 



Controller 2 



I 



M, 



SDF2J A J P^ 



'qSa- + rCtt 
-pSo' + Q'So' 
j-p Ca + a Ca_ 



(Eq. F. 5. 2) 



Controller 3 



^ M SDF3| a J P^ 



A 



- (r + a) So- 
+ (r + a) Ca 
_ pSo- - q Ca_ 



(Eq. F. 5. 3) 



Controller 4 



I 



M SDF4J a J P^ 



( a - r) Scr 
( a - r) Ca 
. p So- + q Cff 



(Eq. F. 5. 4) 



Controller 5 






SDF 5 J A J P^ 



A 



(q - a) Ca 
p Co- - r S a 
_ (q - a) S a _ 



(Eq. F. 5. 5) 



304 



Controller 6 
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(Eq. F. 5. 9) 
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* Note that the equations of this section have been written in 
the form of a column vector because they are simple expres- 
sions with only a few terms. When a number of the above 
equations are added to form a complete control system, the 
form of Equation 3. 2. 3 is recommended. 
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APPENDIX G 



EQUATIONS FOR PARTICULAR MOMENTUM EXCHANGE 
ATTITUDE CONTROL SYSTEMS 

G. 1 INTRODUCTION 

A spacecraft attitude control system of the momentum 
exchange type will usually contain a number of momentum ele- 
ments either of the inertia reaction wheel or of the gyro type. 

In this thesis inertia reaction wheels are referred to as simply 
inertia wheels and have fixed axes of rotation whereas gyros 
are called controllers to differentiate from the conventional 
gyros used in inertial reference systems and in addition to its 
spin motion a controller may have one or more degrees of 
freedom. To determine the moment contributions of a particular 
system consisting of several controllers or wheels as the case 
may be, one needs only to sum the contributions of the individ- 
ual controllers and wheels. A number of the more common 
controllers are represented in Appendix F. 

In this Appendix a number of control systems have been 
defined by choosing various configurations of controllers. The 
following table indicates the controllers which constitute the 
system. 
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List of Control Systems and thei r Cont rollers 
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* The numbers in the parenthesis partially define the control 
system by giving the controllers used for roll, pitch, and 
yaw respectively. Thus the code following the Sun Pointing 
System (0-34-34) indicates that there is no roll control; 
controllers 3 and 4 control both pitch and yaw. A controller 
used in only one axis indicates a single degree-of- freedom 
controller, and for controllers 1 through 6 a controller used 
for two axes indicates a two-degree-of- freedom controllers. 

A controller used for three axes is a three-degree-of- freedom 
controller. 

# These systems are adaptive versions of the (12-34-1234) 
four controller system. 
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Single Controller System (3-3-3) 
Description of System: 



System consists of a number 3 con- 
troller which has three degrees of 
freedom. The rotor of the controller 
is free to precess in two directions 
plus the rotor is capable of being ac- 
celerated about its axis of symmetry. 



Control Logic: 



No control logic is required for the 
single controller system since the 
system contains three degrees of free- 
dom, and the control system input ma- 
trix can be nearly diagonalized by 
suitable arrangement of the matrix. 

+ /3 gives roll to left 
+ y gives pitch nose down 
+ a gives yaw to left 

Moment Contribution of Control System: 
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(Eq. G. 2. 1) 
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S u n Pointing System (0-34-34) 

Description of System: This system is a special case of the 

four controller system. The controllers 
are TDF. 



Control Logic: 
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+ Yo gives pitch nose down 
+ a gives yaw left 



Moment Contribution of Control System: 
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(Eq. G. 2. 2) 
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The equation written in terms of control input variables 
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(Eq. G. 2. 3) 
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Sun Pointing System (12-0-12) 



Description of System: This system is a special case of the 

four controller system. The controllers 
are TDF. 



Control Logic: 
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+ 7^ gives roll right 
+ 0 ^ gives yaw left 



Moment Contribution of Control System: 
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(Eq. G. 2. 4) 
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The equation written in terms of input control variables is 
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(Eq. G. 2. 5) 
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Description of System: Same as four controller system except number 4 controller has failed. 
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Four Controller System (12-34-1234) 



Description of System: 



Four TDF Controllers are mounted 
in positions 1, 2, 3, and 4. 



Control Logic: 
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+ Y^ gives roll to right. 

+ Yg gives pitch nose down. 
+ a gives yaw to left. 
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Moment Contribution of Control System: 
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(Eq.G. 2. 10) 
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The equation written in terms of input control variables is 
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(Eq. G. 2. 11) 
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Six Controller System (56-34-12) 



Description of System: 
Control Logic: 



Six SDF Controllers are mounted 
in positions 1, 2, 3, 4, 5, and 6. 
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+ a^ gives roll to right 
+ ^2 gives pitch nose down 
+ a gives yaw left 

Moment Contribution of Control System: 

-Cov 



M 



(56-34-12)J 



= 2 J p /3 



+ 2J p 3 



0 

0 

0 

-Sa 
+ SQ',- 



1 



0 

-Cov 

0 

+ Sc* 
0 

+ SQ'r 



1 



0 

0 

-Ca'^j 

-Sg^ 

-Sttr 



rj 



y 



(Eq. G. 2. 12) 
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The equation written in terms of control system input variables 
is 
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(Eq. G. 2. 13) 
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Three Controller Orthogonal System (5-3-1) 



Description of System: System consists of SDF Controllers 

Numbers 1, 3, and 5. 



Control Logic: None required since system has 

only three degrees of freedom and 
control system input matrix can be 
nearly diagonalized by arrangement. 

gives roll to right. 

+ 0-2 gives pitch nose down. 

0 

+ gives yaw to left. 



Moment Contribution of Control System: 
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(Eq. G. 2. 14) 
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Controller Zero Momentum System (78-78-1) 
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Two Controller Yaw System (0-0-12) 
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Two Controller Pitch System (0-34-0) 
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(Eq. G. 2. 17) 



Logic 
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(Eq. G. 2. 18) 



Two Controller Roll System (56-0-0) 
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MIT "Vertical Vee” (0-0-9, 10) 



Description of System: System consists of two SDF Controllers. 



Control Logic: 



None proposed since system is to 
be used as passive control system. 
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(Eq. G. 2. 20) 
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Four Controller Sun Pointing System (0-34-12). 



Description of System: Six Controller system without 

Controllers Number 5 and 6. 
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(Eq. G. 2. 21) 
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Three O rthogonal Inertia Wheels (Not a Gyro Syste m) 

Control Logic: None required because system has three 

degrees of freedom and control system 
input matrix can be diagonalized by inspec- 
tion. 
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